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ABSTRACT 
The paper deals with the propagation of its own wave dissipative layered cylindrical bodies interacting with 
deformable (viscoelastic) media. The dynamic behavior of cylindrical bodies is described by equations of 
continuum mechanics. The spectral problem is reduced to solving a system of ordinary differential equations 
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of the first order with variable complex coefficients. Solutions of differential equations are expressed by a 


special cylindrical Bessel functions and Hankel. The frequency equation is solved numerically by Mueller and 
Gauss. The change of own frequency and phase velocity as a function of the wave number. For 
inhomogeneous dissipative mechanical systems found non-monotonic dependence of the imaginary part of 
the phase velocity of the wave numbers. 


Keywords: own waves, dissipative body, Wednesday, spectral problem, the natural frequency and phase 
velocity. 


1. INTRODUCTION 

Regularities of distribution of indignations in continuous environments are of considerable interest 
to many fields of science and technology [1,2]. In this work, unlike known, distribution of own 
waves in layered dissipatively heterogeneous (elements of mechanical system have different 
rheological properties) cylindrical bodies is considered [3,4,5]. 

The analysis of distribution of waves in elastic medium is based on an assumption about justice 
of Hooke's law according to which tension in this point of the environment in timepoints is pro rata 
to deformations in the same timepoints [6,7]. The fact that energy of a series of waves or an 
impulse remains to a constant is a consequence of this assumption. As a rule, influence of these 
"no - gukian" of tension for metals and metal alloys isn't enough, however they are of great 
importance for rubber of similar and construction materials which find essential temporary effects 
under loading [8]. Also different viscoelastic (rheological) properties of dissipatively heterogeneous 


bodies distribution and attenuations of waves significantly influences the law [9,10,11]. 


2. STATEMENT OF THE PROBLEM AND METHODS OF SOLUTION 
In cylindrical system of coordinates { r, 9, z } the distribution of own waves in an isotropic 
viscoelastic body (V) consisting of piecewise and uniform cylindrical bodies of V_k is considered 
(k=1, N), shipped on homogeneous viscoelastic environment (fig. 1). 
Linear the equation of the movtion of mechanical system in a vector form in the absence of 
volume forces takes a form: 
2- 


i,.V7u + (A, +i, )graddivu = 7, a 


,(K =7, 2, 3....N) (1) 


where 
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AF (t)= Ao] FO =|(ReG= oped 


ASO = ba] SO- J Ruta} 2 


f(t)-arbitrary function of time; R,(t—r) ea R,(t—r)-relaxation kernel; A,,,4,9-instant elastic 
moduli(k=1......N); u-displacement vector; p, -density of the medium, k- ordinal number of layers. 
Between the layers is put hard conditions or sliding contact [12] 


r= a, : On = Orr(eHl)? Oa = O-6(«+1)? Ovex = Ore (etl)? 


u, =u B= 9.5 W, = W.y- (3, a) 


KHL? 
If the outer surface is free from multilayer cylinder stresses, whereas 
r=Ayp + O =03 Oyo) =9 Ore =9. (3,6) 


rr(K+1) 
If the displacement vector presented in the form of potential and solenoidal parts 
u=grad p+roty , 


the wave equation have the form 


1 ao 
Ad, -—-— 2 =0; 
A o ot 
1 Oy 
Ay , -+— =0; 
W x eo Ot? 


2, 
ay, Ue Va Te 9, 
‘4 r’- 00 cy, Ot 
Va 20v, 1 OWe 
A + =; 
oe POF 00. 0, OF 
ae 10 0 @ 
= + + af , 
Or? rér 00° 027 
Up = 0g, 4 1 OV 4 OW % : 
or r 00 Oz 
_ 1 , 4 OW x OW y 2 
r00 &~Oor 
u, = 0g, rf OV 4 OW 4 1 OW jeeHL2. al. 
or or or r 06 


Ug, 


where @ — potential of longitudinal waves; Y(W,.Wo> y .) - potential transverse waves. 
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Figure 1 Design scheme of piecewise-homogeneous cylindrical body 


In the construction of representations for motion vector components 


a multilayer cylinder 

VK rea. gehen Ay =e Res oy SF So (& =1,2,...,N)|z,|<00 
we proceed from the equation (4) for the scalar ¢(r,0,z,t)and the vectory(r,0,z,t) potentials. The 
geometry of the object and the natural assumption about the nature of wave motion along the axis 


OZ allow largely to predict the shape of the desired scalar and vector functions. They should 


represent waves traveling along the OZ axis. The solutions of (4) in the form: 


sled life bre 


i —sinn@ 


ia sinn@ ts 
Wik (0.2) Sve Bi fe ” e€ > 
n=0 


—cosn@ 
2 COSNO | sins icy. (5) 
YaltOrzdt)= Sv Bar gh re 


7 sinn@ | sin. jo 
Ve: (-0iz4)= Sv. Je i e€ ) 
n=0 


cosnd 


pase | 9 
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wheren - integer; 7,, - constant wave propagation; @ - natural frequency; pork, z=—. 
ay Ay 
At infinity (r + 0)put Sommerfeld conditions for each component. Substituting (5) into (4) we 


obtain the following ordinary differential equation: 


2 2 


dr? r dr 
avn 1 qa +| B ue m _ 0: (6) 
dr’ be dr Ma =" 


2 
TW da, | 


dr? r ar r 


( NW + 2nNWy Va )BV = 0; 


dy, ldy, 1 
ne = r i . r- (— nv, +2n Wa —Viz Pv. = 0; 


oh 2 2 2 22 oa , 2(1- V ) 
=—t_y?. 20 47) 0.2" gate, 
ie : ue Pr ve . Cox an 1-2v, 


were Q; 
The first two equations in (6) the following solutions for k = 1 of the first (solid) and the outer 
cylinder (k=N): 


FY, (a,r), k a 1; 
= 7, 
?, (r) cae = N: (7,a) 


M,,J, (B,r ), k =1; 
y(r)= (7,6) 


Mi,H(B,r), k=N; 
where J, (@,) - Bessel functionn - th order. H'(Z,r) - Hankel function of the first kindn - th order. 
The decision of the other two equations in (6) is also expressed in terms of Bessel functions 

and Hankel: 


ok 1 AP Land nt k 
v(r)= | ie ae ee 


Diggit n- O (Br )+ UH) (2.7), k= N; 
Solutions (7a), (7b) and (8) of the system of differential equations (6) will contain a 6k-2 
arbitrary constants. | take this opportunity to go to a large extent arbitrary values when choosing a 


permanent believe further L,, = LZ, =0,1.¢.v, =Wq- 


Moving k cylinder is expressed in terms of Bessel and Neumann functions of n-th order 
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U,. = > {rl 1a (7,7) + Arlee )]+ [Ast (72,17) + ve & (7.,.7)|- 


(9) 


= a Aaa (7.,r)+ Auli cosngel") 


2K 


Ua = = Vx Avent, (Vict) + AdenY, (v,.7)|+ _ [Ay,,.J7 (Yar) + Aan¥(%ou0)|- 
n=0 n 


-— [Asiad (Yo, r) + Age Ny (Yo. nl sin nge re) 


2k 


0 ) —_— 
Ux = » in Ain! (7yr)+ Ay,Y, (rar TAs (7,7) + Asin’, (7,,r)eosnge" mae ’ 


n=0 


where 
a _ 7) = 
ee aes ae rk Vi. = fae = 
Hy, = Oley; cf = y=malil, (m =1,2,...) (k =1,2,3) 


To determine the arbitrary constants A,,,,A5jn> Asin> Adin? Asin? Aon USEd the boundary conditions 


(3) .Then we obtain a system of homogeneous algebraic equations6x +3 unknown and6x +3 
equations. A necessary and sufficient condition for the existence of solutions of this system, is the 
vanishing of its determinant. The order of the primary identifier (6x +3x6x+3) and elements are 
expressed in terms of Bessel and Neumann functions of n-th order complex argument. This 
equation gives the dispersion equation for the dissipative systems. 


Under the dispersion characteristics are understood according to the phase (C =C, +iC,) 
and group velocityV =V, +iV, the wave number (y,,) while various parameters of the mechanical 


system. It is known that the magnitudes of C and V related to the value of the root of the dispersion 
equation 


A(C,@,7,,,4) =0, (10) 
where w- complex frequency; y,-wave number, 7 — wavelength. The phase and group velocity is 


related to the root meaning of the dispersion equation some complex dependencies. Thus, to be 
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able to calculate the dispersion characteristics, it is necessary to make a study of the roots of the 
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equation (10) at the points of the complex plane, and also to develop a method of numerical 


determination. 

The work for the solution of the transcendental equation (10) applies the method of Mueller, at 
each iteration of the method applied by Muller Gaussian with the release of the main element. 
Thus, the solution of equation (10) does not require disclosure of the determinant. As an initial 
approximation we choose the phase velocity of the waves corresponding elastic system. For 


complex roots Muller method simplifies the calculations and provides faster convergence [13,14]. 


3. THE PROPAGATION OF TRANSVERSE WAVES IN AN INFINITELY LONG CYLINDRICAL 
SHELL, LOCATED IN AN ELASTIC MEDIUM 

The main objective of the research study of the existence of the phase velocity of 
propagation of the geometric and physical and mechanical parameters of the system. The basic 
equations of the theory of elasticity for such tasks are reduced to a plane problem (U, =U, = 0). 


The components of the displacement vector in the cylinder and its environment are represented as: 


tn = alfa! LHe Knee! Lae K|| e : 
Or dr dr : 


U 


Ugg =—C" fk (Kne* ; 
dr 


dr dr — 
Y re =— 1 F F e 
fal HO (K,r)+B' ate ia) 
r 


dr 


I d° Ore I d° Qe ~K;z 
Ay ne —— a A ie + 
-iK 


To determine the arbitrary constantsA!, B'andC'' It used the boundary conditions r = ai: ore = 0 
and r = a2Uei= Ue2z, ORo1 = GO Re2We Obtain a homogeneous system of algebraic equations of the 
third order. Their conditions of existence of nontrivial solutions we obtain the following dispersion 


equation: 
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-K?H®(K,a,) -K, H®(K,a,) 0 
-K, HO(K, a,) ~K, H(K,a,) ~ KK, (K,,) =0, (11) 
—K? MAO (K, a») — KP MHOC, a,) -K,u,K, (K,a,) 
2 2 
K,= a — Kz > Ky — = »Ky = a ? Co ~ oe = £2 
s1 s2 Pi P2 
when K> >K*1. 
V 
U [m /cek] 
Cp2 ” 
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U 
80 
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Figure 2 Changing the phase and group velocities, depending on the wavelength. 


Cg1 = 200 m/s; Cg2= 1100 m/s 
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V U [m /cek] 
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Cp1 
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Figure 3 Changing the phase and group velocity as a function wavelength from 


The numerical results are shown in Figure 2. and Figure 3. Note that an increase in the thickness of 
the layer of the first and second mode phase velocity gradually decreases. From the figures it can 


be seen that the phase and group velocity in the area of wavelengths in the range of C,,(A(C,,,. In 


the of short waves phase velocity may be less than the velocity of shear waves. 


Distribution of natural waves in the three-layered cylindrical body (Figure 1). Internal 
(r=a7) and external (r=a4) the surface of the three-layer cylinder is free from stress. Then the order 
of the determinant of the system of homogeneous algebraic equations (18x18), which is the 


dispersion equations (a, —a, =Ah,, a,-a, = Ah,, a,-a, =Ah,). 
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Figure 4a Changing the complex natural frequencies of the wave number (Dissipative 


homogeneous mechanical system). 


In the calculations we take the following values: 


Ah; = 0,02 m, Ah2 = 0,03 m, Ah3 = 0,02 m, E, =5,5.10°Pa, E, =3,5.10° Ma, E; =5,5.10°Pa p, =27 


Kg/m>. p, =11kg/m*. p, =27kg/m?. A=0,048; £B=0,05; a=0,1. 


Baez 5 


ARTICLE 


Consider two options for a dissipative system. In a first embodiment, dissipative system is 
structurally homogeneous. 


t.. 


Figure 4b Changing the complex natural frequencies of the wave number (Dissipative 


heterogeneous mechanical system). 
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The dimensionless wave number it varies between 0 - 2. The calculation results are shown in 


Figure 4 as well. The dependence of the real and imaginary parts of the phase velocity of the 


dimensionless wave numbery,, It was monotonous, and dependence is the same for the real and 


imaginary part of the phase velocity. In the second embodiment, a dissipative system is structurally 
inhomogeneous, i.e. Rheological properties of the core layer (n = 2) is equal to zero, the other 
parameters the same as adopted above. The calculation results are shown in Figure 4b. The 


dependence of the real part of the phase velocity of the dimensionless wave number of 7, It is the 


same as for the homogeneous system: corresponding curves coincide up to 5%. Dependence of 


the imaginary part of the phase velocity of they, nonmonotonic it turned out. 


Of particular interest is the minimum valueéa fixed damping coefficient: 


Sq@=min (- M %), Sc=min (-cx) k= 1,2,.....K,(10) 


here 5: — coefficient determining the damping properties of the system. 


For a homogeneous system coefficient 5, entirely determined by the imaginary part of the first 
modulo complex phase velocity. For inhomogeneous systems as a factor 6< can act as an imaginary 
part of the first and the second frequency depending on their values. "Turn the Tables" occurs 


when the characteristic value y,,, while the value of the real parts of the first and second frequencies 


are closest. 


Coefficient 6 @ and 5, at the specified characteristic value has a pronounced maximum. 


4. CONCLUSIONS 
1. In the non-monotonic dependence of the imaginary part of the solution to the problem of wave 
propagation in inhomogeneous media dissipativno- discovered natural frequencies (or imaginary 
part of the phase velocity) of physical - mechanical and geometric parameters of the system. 

2. It was found that the phase velocity of the highest forms of extension and torsion waves exceed 


the highest possible speed with the waves in an infinite medium, the group velocity never exceeds 
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S.Takzheuotanavili that nondispersive medium grupovogo rate of 10-15% previschaet comparison 


dispersion medium. In other words, the shape of the pulses as they propagate not remain the same 


as in homogeneous elastic bodies. 


3. Increase in viscosity reduces the real and imaginary parts of the phase velocity (or real and 


imaginary parts of the natural frequencies) of up to 20% for the homogeneous dissipative 


mechanical systems. For inhomogeneous dissipative mechanical systems increase in viscosity 


reduces the real part of the phase velocity (or the real part of the natural frequencies) to about 


18%, and their imaginary parts, changed radically. 
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